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We discuss the dependence of observables on the chemical potential in 't Hooft's large- QCD. 
To this end, we use the worldline formalism to expand the fermionic determinant in powers of 1/A'^. 
We consider the hadronic as well as the deconfining phase of the theory. We discuss the origin of the 
sign problem in the worldline approach and elaborate on the planar equivalence between QCD with 
a baryon chemical potential and QCD with an isospin chemical potential. We show that for C-even 
observables the sign problem occurs at a subleading order in the expansion of the fermionic 
^ , determinant. Finally, we comment on the finite A'' theory. 
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I. INTRODUCTION 
• 

Understanding the QCD phase diagram is an important and longstanding problem (see Ref. [l| for a review) . Since 
' the theory at a generic point in the (T, /i) plane is strongly coupled, one has to resort to lattice simulations. It is, 
however, difficult to carry out Monte Carlo simulations in a theory with a chemical potential due to the so called 
■ "sign problem" . When a theory admits a nonzero chemical potential, the fermionic determinant becomes generally 
] complex, and large oscillations of the phase make the Monte Carlo simulation practically impossible Q • 
, We could expect that the sign problem simplifies in the 't Hooft large- iV limit of QCD (Sf. In this limit, fundamental 
matter is quenched, and the hope is that the problem which arises from the determinant is a subleading l/N effect. 
, The problem was first discussed at the diagrammatic level by Cohen in Ref. A related idea which was proposed 
recently is the large- equivalence between QCD and other theories which admit a real fermionic determinant [5|-|8|. 
' In this paper we use the worldline formalism [o^ to study the dependence of observables on the chemical potential 
' in large-A^ QCD. In this approach, the fermionic determinant is written as a sum of Wilson (or Polyakov) loops. 
' Moreover, the expansion in loops corresponds to a systematic expansion of the fermionic determinant in powers of 
l/N^. The main disadvantage of this approach is that it is perturbative in l/N, and hence the expansion is around 
the Yang-Mills vacuum (without quarks). As a result, the discussion will be restricted to the phases of the theory 
with small fj, where there is no breaking of baryon (or isospin) number. 

The main goal of the paper is to understand how observables depend on the chemical potential in the large- A'^ theory. 
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, We will also discuss the "sign problem" . We will explicitly decompose the determinant into real and imaginary parts 
and shed light on the origin of the problem. The outcome of our analysis is that while the large- TV theory admits 
a sign problem, the problem is not as severe as in the finite-A'^ case. More precisely, we will show that to leading 
fermionic contribution in N, C-even (real) operators do not couple to the phase of the fermionic determinant. In 
other words, for C-even operators, the sign problem is a subleading effect. 

This understanding will be then used to clarify the meaning of the recently proposed large- A^ equivalence between 
QCD theory with a baryon chemical potential and a QCD theory with an isospin chemical potential [1, Q. 

We conclude the paper with a discussion about the implications of our analysis on the finite- A^ theory. 
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II. THE 1/N EXPANSION VIA THE WORLDLINE FORMALISM 

Let us review the worldhne formahsm • Consider SU (N) gauge theory with one flavor and without a chemical 

potential. We formulate the theory on R'^ x S^, where is a temporal circle with a radius R = 1/2ttT. The basic 
idea is to express the fermionic determinant in terms of Wilson loops. The precise relation between the Euclidean 
fermionic determinant and Wilson loops is as follows: 



and 



where 



det(^ + m) =expr[A], (1) 



Z = J VA^eM~SYM)eMnA]) , (2) 



X J VxVijj (-^ly exp |- dr Qi^i^ + ^^A^^^'j | 



X TrPexpjz dr (^A^i^ - -V^F^^^^j T j , (3) 

with a;^(0) = x^{t) and uj being the winding number of the path around the compact direction, {t"} is a set of 
matrices which transform in the fundamental representation of the SU{N) Lie algebra. The factor (—1)" implements 
the antiperiodic boundary conditions for fermions in the worldline formalism. Thus F, is a sum over (super-)Wilson 
loops. The sum is over contours of all sizes and shapes. 
Let us introduce a chemical potential /i: 

det {]/> + m) ^ det {p + m + . (4) 

In the worldline formalism, a chemical potential is easily implemented by introducing a constant background U{1) 
gauge connection 

{A,) = -^5"^^,, (5) 

resulting in the introduction of a factor 

exp{fi J dx°) = exp{j;uj) (6) 

since the integral along the a;° direction (the circle) measures the winding uj. Let us decompose T accordingly to the 
various topological sectors which correspond to the winding number 

oo oo 

r(M)= r"^HM)= E r(")(M = 0)exp(^u;). (7) 



a; — — oo a; — — oo 



While we cannot prove that the above expansion ([7]) is always convergent, it is expected to converge when m/fi is 
large enough, since the quark mass serves as an IR regulator such that each closed loop carries a factor ^ exp(— ^ |w|). 
The free energy at generic N is defined as 

F=-T log Z = Fym - T log(er)YM ■ (8) 

We wish to consider now the expansion of expF. As was noted in Ref. [Tlj, the connected fc-point function of F is 
suppressed by powers of N: 

k times 
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Therefore, the expansion of expF in powers of T is an expansion in l/N. The approximation expF ^ 1 corresponds 
to quenching. In order to consider the effect of fermionic matter, one has to go beyond the quenched approximation, 
namely, to consider one Wilson loop (which corresponds to one power of F) 

expF-l + F. (10) 

In particular, the leading 0{N) contribution of ferniions to the free energy is 

oo 

F/ = -Tlog(e^)YM--T(F)YM = -r ^ exp(^aj)(F('^)(/i = 0))ym • (11) 

UJ— — 00 

Note that the expectation value is calculated in the pure Yang-Mills vacuum. Next, we consider a generic Wilson 
loop O^^ ^ with winding number uj' . At 0{N)^ the expectation values of O^^ ^ are given by 

oo 

(0("'))qcd = (0^"'))ym+ exp(^..)(O(-')F(-)(/x = 0))YM. (12) 

cj— — oo 

In the next sections, we will elaborate on the implications of the expression (I12p . 



A. 11 dependence in the confining phase 

Let us restrict our attention to the hadronic confining phase. In the confining phase, the expectation value of any 
loop wrapping uj times around the circle vanishes unless lj mod = 0. At the leading fermionic contribution in N , we 
can neglect 0{N) windings. Therefore, the free energy PT|) will get a contribution only from the sector with w = 0, 

F/--T(F("=")(m = 0))ym. (13) 

Thus, we conclude that in the hadronic phase the leading fermionic large-iV contribution to the quantity F/T is /i 
and temperature independent, 



(14) 



T=0,/i=0 

The same conclusion holds for any Wilson loop with a zero winding number, 

(O("'=o))(r^^) ^ (O("-"))(r = 0,Ai = 0). (15) 
Let us consider a Wilson loop O*^" with lo' ^ 0. In this case Eq. (|12p is saturated by w' + w = 0, namely, 

(0^"'')qcd = (0("'))ym + cxp(-^c.')(0("')r(-"''(A* = 0))ym . (16) 

Even in the confining phase, expectation values of Wilson loops with a nonzero winding admit a nontrivial dependence 
on the chemical potential. 



B. ^ dependence in the deconflning phase 

In the deconfining phase, the Polyakov loop expectation value is nonvanishing. Therefore the free energy is not 
saturated by a Wilson loop with w = 0, but it rather acquires contributions from all the topological sectors. In 
particular, the free energy is ii and temperature dependent. Using Eq. ([TT|) and charge-conjugation invariance of the 
Yang-Mills vacuum, we arrive at 

Ff = -r(F(o)(/i = 0))ym - 2T^(F(")(/. - 0)}YMCOsh(^c^) . (17) 
In general, observables in the deconfining phase are both and temperature dependent. 
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C. A comment on mesonic correlators 

An important question is how the mesonic spectrum is affected by the chemical potentiaL For simphcity, we will 
consider the case of a theory with a single flavor. Our discussion applies also to a multiflavor theory with a baryon 
chemical potential. 

The mesonic spectrum of a large- iV theory is extracted from a two-point function as follows 

(^vE'(x) , ^^{y)) = J ^ exp(*fc(x - y))G{k') (18) 



with 



A:2 + M2 ' 

where /„ are coupling constants and Af„ are the meson masses. In order to compute the two-point function (|18p in 
the worldline formalism, one has to add to the action a source / d'^x J{x)^'i>{x) and to differentiate the partition 
function twice with respect to the current J. The outcome [T^, to leading order in N, is that Eq. is given by a 
summation over all Wilson loops which pass via x and y, 

{^^{x),^^{y))^{^,^y{^,)). (20) 

Each Wilson loop in the sum carries a factor of exp(/ia;/T) where w is its winding number. 

However, in the confining phase, the correlator is saturated by Wilson loops with a zero winding number: hence, 
we conclude that the leading large- meson spectrum is not affected by the presence of the chemical potential. 

III. THE ORIGIN OF THE SIGN PROBLEM 

In this section, we elaborate on the origin of the sign problem in the worldline approach. 
Let us start by proving the following lemma, valid at fixed gauge configuration: 

r('^)(^ = o)* = r(-"^'(/x = o) (21) 

Proof: consider the theory with a pure imaginary chemical potential. In this case, 

oo oo 

J2 r('^)(*M)= E r(-)(0)exp(i^L.)=logdet(^ + m-fiAi70) , (22) 



T 

LU — — 00 LJ — — 00 



hence, 



r(-)(0) = ['^^ ^ exp(-z^u;) logdet + m + ^^^7°) • (23) 







Since the determinant is positive (hence, the logarithm is real), taking the complex conjugate is equivalent to replacing 
Lo — > —Lo. Equation ([23]) can be written also in the following form: 

rM(0) + r(-")(0) = 2 f -^cos(^w)logdet(^-Hm-f i^7°) (24) 
J Zirl 1 

r(-)(0)-r(-)(0) = -2z j ^sin(|u;)logdet(^ + m-f *M70) , (25) 

which separates the real and imaginary parts. In particular, note that r(°)(0) is real. 

The above analysis has implications on the origin of the sign problem. When /i = 0, the fermionic determinant is 
real because we sum over all configuration with all winding 

oo 

det (^ + m) = exp ^ V^'^^ (m = 0) , (26) 
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and in spite of r^'^)(/x = 0) being complex, the imaginary part cancels between F^"^) and F^""^. However, when we 
introduce a real chemical potential, F^"^ is weighted by exp(^a;), while F'^"'^) is weighted by exp(— ^w), and the 
imaginary part no longer cancels: 

oo 

^ F(-)(a.) = r(o)(0) + J2 (exp(^^)F(")(0) + exp(-^c.)F(-)(0)) (27) 

uj—~oo uj>0 

This is the reason for the sign problem. From Eq. , we learn that 

5R(F(")) = 5R(F(-")) (28) 
3(f(")) = -S(f(-'^)), (29) 

and the fermionic determinant becomes 

det {p + m + Ai7°) = exp <j F^*') (0) + 2 ^ cosh(^w)5R(F(") (0)) I x 



i<;>0 



exp|2z5]sinh(^c.)S(F(")(0))| . 



(30) 



In particular, note that the notorious phase which causes the sign problem is given by 

2^sinh(^..)S(F(-)(0)). (31) 

In the large- limit, in expectation values of observables with a well-defined large- iV limit, the fermionic determinant 
can be replaced by 



det {p + m + i^Lj°) 1 + r(°) (0) + 

(^..)3?(F(-)(0)) + 2zs-_^ 



J2 (2cosh(^..)3?(F(-)(0)) + 2zsinh(^u;)3(F(-)(0))) . (32) 



cc;>0 



We want to investigate now for which observables the imaginary part in the previous equation does not contribute 
and can be dropped. 

The free energy in the confining (hadronic) phase is saturated by F(°'(/i), which is real, so obviously there is no sign 
problem in this case: only configurations with zero winding number contribute to the free energy. A more interesting 
case is the deconfining phase. The imaginary part in Eq. ((32|) drops once we insert it into Eq. pT)) . because the 
Yang-Mills vacuum is invariant under charge-conjugation ((F'^'^^ (0))ym is real). 

Let us consider now observables which are Wilson loops of arbitrary shape with generic winding number w'. At the 
first leading fermionic contribution: 

oo 

(0("')) = (0(-'))ym + 2 J2 exp(|c.)(O(-')F(")(^ = 0))c,YM = 

UJ— — 00 

oo 

= (0("'))ym + 2 cosh(^w)(KO("') KF(")(/X = 0))e.YM + 

CJ— — OO 

oo 

-2 J2 sinh(^c^)(30("') 3F(")(M = 0))e,YM. (33) 

In the confining phase and for Wilson loops with no winding oj' = 0, only T^^^{ii), which is real, contributes because 
of center symmetry. In the deconfined phase, the imaginary part of the determinant generally contributes. 

However it is clear that the imaginary part of the determinant couples only to the imaginary (or C-odd) part of 
O^'^ \ If one consider only the C-even part, then the imaginary part of the determinant decouples from any generic 
connected expectation value (in both phases), thanks to the C-invariance of the YM vacuum: 

{m[^'''^ • ■ • Koi""^)c = (3?o["^^ • • • Koi""^)c,YM + 

oo 

2 ^ cosh(^w)(»oi"^^---3?oi""^ 5RF(")(^ = 0))c,YM. (34) 
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The phase of the determinant decouples from 5R0*^" ^ at the leading fermionic contribution, which is the first subleading 
contribution in The sign problem is therefore a l/N"^ effect for 3?0^" 

In the confining phase only, the imaginary part of the determinant decouples also from the loop 0^°-' with zero 
winding number, since only the sector with zero winding contributes (because of center symmetry), and this is /i- 
independent: 

(0(°) • • • OW)e(0^°^ • • • 0("))e,YM + 2(0^°) • • • 0^^)V(^\^, = 0))e,YM ■ (35) 

Summarizing, the sign problem does not exist at the leading fermionic contribution in for a class of observables 
which include real parts of Wilson loops with any winding in both phases, and Wilson loops with no winding in the 
confining phase only. 



IV. LARGE-A/^ EQUIVALENCE BETWEEN BARYON AND ISOSPIN CHEMICAL POTENTIALS 

Let us consider an SU{N) gauge theory with two identical massive flavors. Let us assign the same chemical potential 
fi to both flavors ("baryon chemical potential"). After an integration of the Fermi fields, we obtain a square of the 
single- flavor fermionic determinant, 

(dct(^ + 771 + ^7°))' • (36) 

In the worldline formalism, implementing two identical flavors is achieved, with respect to the single flavor, by the 
assignment 

exp(r(Ai)) exp(2r(M)) . (37) 

Alternatively, one can assign a chemical potential fi to one flavor and a chemical potential —fi to the other flavor 
("isospin chemical potential"). In this case, we obtain a product of fermionic determinants, 

(det(^ + m + /i7°)) (dct(^ + m - ^7°)) , (38) 

resulting in 

exp(r(/i)) ^ exp(r(Ai) + r(-/i)) . (39) 

It has recently been argued [^, Q that the above two theories are planar-equivalent in the hadronic and deconfining 
phases. This means that a common sector of observables exits in the two theories, such that all the connected 
expectation values of such observables are identical in the large- TV limit in the two theories. An earlier perturbative 
argument in favor of the equivalence was given in Ref. jj]. The equivalence and its breaking at smaller N are also 
supported by a recent numerical lattice simulation [l^ . 

Considering an isospin chemical potential instead of a baryon chemical potential amounts to dropping the imaginary 
part of the determinant p2p at the leading nontrivial order in 1/N. Indeed, the fermionic determinant of a two- flavor 
theory with a baryon chemical potential is given by 

(det(.^ + ))'^l + 2r(°)(0)-f 

+2^ (2cosh(^^)3?(r(")(0)) + 2isinh(^w)3(r(")(0))) , (40) 

while the fermionic determinant of a two-flavor theory with an isospin chemical potential is given by 

det(^ + m + /i7°) det(^ + m - ^7°) ^ 1 + 2r(°)(0) -f 

+2 5^(2cosh(|a;)K(r(")(0))) . (41) 

Therefore, all the observables analyzed in the previous section which decouple from the imaginary part of the 
determinant with baryonic chemical potential belong to the common sector. Put differently, the common sector of 
the two theories includes gluonic C-even operators. For these operators, the sign problem is a subleading effect - they 
do not couple to the phase of the fermionic determinant (which exists only for the theory with a baryon chemical 
potential) at the leading fermionic contribution. In addition to these C-even operators, there exist other quantities 
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which belong to the common sector: for example, the charge density which is obtained by differentiating the free 
energy with respect to the chemical potential. Obviously, derivatives with respect to /i of quantities which belong to 
the common sector also belong to the common sector. 

The outcome of the above analysis is that in the large- limit of QCD, in a limited part of the phase diagram and 
for certain operators, it is justified to replace the theory with a baryon chemical potential by a theory with an isospin 
chemical potential. 



V. DISCUSSION 



In this paper we focused on the large- limit of a QCD theory with a chemical potential. 

Let us discuss what happens when we introduce higher 1/iV corrections. The simplest case is the free energy in the 
hadronic phase, 

F/T = ^Vm/T - log ( V i- V exp(^Z7V)(r("^)r("^)...r("^)) ) (42) 



Coo ^ oo 
fc=0 i=-oo 



with + UJ2 + ■■■ + oJk = IN, where / is an integer. Namely, for a given k, one has to consider all the partitions of N, 
2N, 3N and so on. These configurations contribute to the free energy and lead to a nontrivial /i dependence [13, E^- 
For this reason, our conclusions about the large- iV theory cannot apply for the finite- theory. 

Another important issue is the sign problem at finite A''. At infinite A^, for C-even observables, only the real part 
of r contributes: hence, expF is positive, and the measure is positive-definite. This is not the case at finite N and 
the sign problem is severe. 

The main result of our paper is a field-theory understanding of the planar equivalence between the theory with a 
baryon chemical potential and the theory with an isospin chemical potential. It is now clear why it is useful to carry 
out lattice simulations using an isospin chemical potential: there is a limit of QCD (the large- A^ limit) where this 
procedure is justified. At finite A^, and, in particular, for SU(3), we expect an error of 1/N in general, and an error 
of 1/A^^ for purely gluonic C-even observables. 

Among future directions of investigation we would like to mention planar equivalence with SO{N) theories and 
theories with adjoint matter. It is straightforward to generalize our discussion to these cases. 
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